Mac Power 1997 January

home *** CD-ROM | disk | FTP | other *** search
/ Mac Power 1997 January / macpower199701.bin / AMUG / Programming_10 / Videotoolbox 96-08-22-c.sit / VideoToolbox / (Utilities) / Quick3 / TestPsychometricFit.c < prev next >
MacBinary | 1995-06-17 | 5.0 KB | [TEXT/CWIE]
open in:
MacOS 8.1 extracted |
Win98 extracted |
DOS extracted
view JSON data |
view as text

This file was processed as: MacBinary (archive/macBinary).

You can browse this item here: TestPsychometricFit.c
Confidence	Program	Detection	Match Type	Support
`66%`	dexvert	Compact Compressed (Unix) `(archive/compact)`	ext	Supported
`10%`	dexvert	MacBinary `(archive/macBinary)`	fallback	Supported
`1%`	dexvert	Text File `(text/txt)`	fallback	Supported
`100%`	file	MacBinary II, inited, Sat Jun 17 01:31:23 1995, modified Sat Jun 17 01:31:23 1995, creator 'CWIE', type ASCII, 4389 bytes "TestPsychometricFit.c" , at 0x11a5 410 bytes resource	default (weak)
`99%`	file	data	default
`74%`	TrID	Macintosh plain text (MacBinary)	default
`25%`	TrID	MacBinary 2	default (weak)
`100%`	siegfried	fmt/1762 MacBinary (II)	default
`100%`	lsar	MacBinary	default
id metadata
key	value
macFileType	`[TEXT]`
macFileCreator	`[CWIE]`
hex view
+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 15 54 65 73 74 50 73 | 79 63 68 6f 6d 65 74 72 |..TestPs|ychometr|
|00000010| 69 63 46 69 74 2e 63 00 | 00 00 00 00 00 00 00 00 |icFit.c.|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 54 45 58 54 43 57 49 | 45 01 00 00 00 00 00 00 |.TEXTCWI|E.......|
|00000050| 00 00 00 00 00 11 25 00 | 00 01 9a ac 08 19 2b ac |......%.|......+.|
|00000060| 08 19 2b 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |..+.....|........|
|00000070| 00 00 00 00 00 00 00 00 | 00 00 81 81 71 b6 00 00 |........|....q...|
|00000080| 2f 2a 0d 54 65 73 74 50 | 73 79 63 68 6f 6d 65 74 |/*.TestP|sychomet|
|00000090| 72 69 63 46 69 74 2e 63 | 0d 43 6f 70 79 72 69 67 |ricFit.c|.Copyrig|
|000000a0| 68 74 20 28 63 29 20 31 | 39 39 30 2d 31 39 39 32 |ht (c) 1|990-1992|
|000000b0| 20 44 65 6e 69 73 20 47 | 2e 20 50 65 6c 6c 69 0d | Denis G|. Pelli.|
|000000c0| 54 68 69 73 20 69 73 20 | 61 20 73 69 6d 70 6c 65 |This is |a simple|
|000000d0| 20 64 72 69 76 65 72 20 | 74 6f 20 73 68 6f 77 20 | driver |to show |
|000000e0| 74 68 61 74 20 50 73 79 | 63 68 6f 6d 65 74 72 69 |that Psy|chometri|
|000000f0| 63 46 69 74 2e 63 20 77 | 6f 72 6b 73 2e 0d 0d 49 |cFit.c w|orks...I|
|00000100| 20 75 73 65 20 74 68 65 | 20 57 65 69 62 75 6c 6c | use the| Weibull|
|00000110| 20 70 73 79 63 68 6f 6d | 65 74 72 69 63 20 66 75 | psychom|etric fu|
|00000120| 6e 63 74 69 6f 6e 20 74 | 6f 20 70 72 6f 76 69 64 |nction t|o provid|
|00000130| 65 20 74 68 65 20 70 72 | 6f 62 61 62 69 6c 69 74 |e the pr|obabilit|
|00000140| 79 20 61 74 20 65 61 63 | 68 20 63 6f 6e 74 72 61 |y at eac|h contra|
|00000150| 73 74 0d 61 6e 64 20 74 | 68 65 20 4e 75 6d 65 72 |st.and t|he Numer|
|00000160| 69 63 61 6c 20 52 65 63 | 69 70 65 73 20 62 6e 6c |ical Rec|ipes bnl|
|00000170| 64 65 76 28 29 20 62 69 | 6e 6f 6d 69 61 6c 20 64 |dev() bi|nomial d|
|00000180| 65 76 69 61 74 65 73 20 | 66 75 6e 63 74 69 6f 6e |eviates |function|
|00000190| 20 74 6f 20 73 69 6d 75 | 6c 61 74 65 20 74 68 65 | to simu|late the|
|000001a0| 20 61 70 70 72 6f 70 72 | 69 61 74 65 0d 6e 75 6d | appropr|iate.num|
|000001b0| 62 65 72 20 6f 66 20 74 | 72 69 61 6c 73 20 61 74 |ber of t|rials at|
|000001c0| 20 74 68 61 74 20 63 6f | 6e 74 72 61 73 74 2e 0d | that co|ntrast..|
|000001d0| 0d 54 68 65 6e 20 49 20 | 63 61 6c 6c 20 50 73 79 |.Then I |call Psy|
|000001e0| 63 68 6f 6d 65 74 72 69 | 63 46 69 74 28 29 20 61 |chometri|cFit() a|
|000001f0| 6e 64 20 61 73 6b 20 69 | 74 20 74 6f 20 66 69 74 |nd ask i|t to fit|
|00000200| 20 74 68 65 20 57 65 69 | 62 75 6c 6c 28 29 20 66 | the Wei|bull() f|
|00000210| 75 6e 63 74 69 6f 6e 20 | 74 6f 20 74 68 65 20 64 |unction |to the d|
|00000220| 61 74 61 2c 20 61 6e 64 | 0d 74 68 65 6e 20 49 20 |ata, and|.then I |
|00000230| 70 72 69 6e 74 20 6f 75 | 74 20 74 68 65 20 72 65 |print ou|t the re|
|00000240| 73 75 6c 74 73 2e 20 41 | 73 20 6f 6e 65 20 77 6f |sults. A|s one wo|
|00000250| 75 6c 64 20 65 78 70 65 | 63 74 2c 20 74 68 65 20 |uld expe|ct, the |
|00000260| 66 69 74 73 20 61 72 65 | 20 76 65 72 79 20 67 6f |fits are| very go|
|00000270| 6f 64 2e 20 49 74 20 73 | 65 65 6d 73 0d 72 65 61 |od. It s|eems.rea|
|00000280| 73 6f 6e 61 62 6c 65 20 | 74 6f 20 72 65 6a 65 63 |sonable |to rejec|
|00000290| 74 20 66 69 74 73 20 61 | 74 20 74 68 65 20 35 25 |t fits a|t the 5%|
|000002a0| 20 73 69 67 6e 69 66 69 | 63 61 6e 63 65 20 6c 65 | signifi|cance le|
|000002b0| 76 65 6c 2e 20 53 6f 20 | 65 78 70 65 63 74 20 74 |vel. So |expect t|
|000002c0| 6f 20 72 65 6a 65 63 74 | 20 61 62 6f 75 74 0d 35 |o reject| about.5|
|000002d0| 25 20 6f 66 20 79 6f 75 | 72 20 66 69 74 73 20 65 |% of you|r fits e|
|000002e0| 76 65 6e 20 69 66 20 61 | 6c 6c 20 74 68 65 20 61 |ven if a|ll the a|
|000002f0| 73 73 75 6d 70 74 69 6f | 6e 73 20 6f 66 20 74 68 |ssumptio|ns of th|
|00000300| 69 73 20 6d 6f 64 65 6c | 20 61 72 65 20 63 6f 72 |is model| are cor|
|00000310| 72 65 63 74 2e 0d 0d 4f | 70 74 69 6f 6e 61 6c 6c |rect...O|ptionall|
|00000320| 79 2c 20 61 6c 73 6f 20 | 63 6f 6d 70 61 72 65 73 |y, also |compares|
|00000330| 20 51 55 49 43 4b 33 20 | 77 69 74 68 20 51 55 45 | QUICK3 |with QUE|
|00000340| 53 54 2e 20 28 4e 6f 74 | 65 3a 20 51 55 45 53 54 |ST. (Not|e: QUEST|
|00000350| 20 69 73 20 6e 6f 74 20 | 70 61 72 74 20 6f 66 20 | is not |part of |
|00000360| 74 68 65 20 0d 56 69 64 | 65 6f 54 6f 6f 6c 62 6f |the .Vid|eoToolbo|
|00000370| 78 2e 20 53 6f 72 72 79 | 2e 29 0d 0d 48 49 53 54 |x. Sorry|.)..HIST|
|00000380| 4f 52 59 3a 0d 34 2f 36 | 2f 39 30 09 64 67 70 09 |ORY:.4/6|/90.dgp.|
|00000390| 77 72 6f 74 65 20 69 74 | 2e 20 53 65 65 6d 73 20 |wrote it|. Seems |
|000003a0| 74 6f 20 77 6f 72 6b 20 | 66 69 6e 65 20 66 6f 72 |to work |fine for|
|000003b0| 20 61 6c 6c 20 63 61 73 | 65 73 2c 20 31 20 74 6f | all cas|es, 1 to|
|000003c0| 20 31 30 30 2c 30 30 30 | 2c 30 30 30 20 74 72 69 | 100,000|,000 tri|
|000003d0| 61 6c 73 20 70 65 72 0d | 09 09 09 63 6f 6e 74 72 |als per.|...contr|
|000003e0| 61 73 74 2c 20 61 6e 64 | 20 66 65 77 20 61 6e 64 |ast, and| few and|
|000003f0| 20 6d 61 6e 79 20 63 6f | 6e 74 72 61 73 74 73 2e | many co|ntrasts.|
|00000400| 0d 31 31 2f 31 38 2f 39 | 32 20 64 67 70 20 61 64 |.11/18/9|2 dgp ad|
|00000410| 64 65 64 20 63 6f 6d 70 | 61 72 69 73 6f 6e 20 77 |ded comp|arison w|
|00000420| 69 74 68 20 51 55 45 53 | 54 2e 0d 32 2f 32 30 2f |ith QUES|T..2/20/|
|00000430| 39 33 09 64 67 70 09 61 | 64 64 65 64 20 63 61 6c |93.dgp.a|dded cal|
|00000440| 6c 20 74 6f 20 52 65 71 | 75 69 72 65 28 29 2e 0d |l to Req|uire()..|
|00000450| 39 2f 35 2f 39 34 20 64 | 67 70 20 72 65 6d 6f 76 |9/5/94 d|gp remov|
|00000460| 65 64 20 61 73 73 75 6d | 70 74 69 6f 6e 20 69 6e |ed assum|ption in|
|00000470| 20 70 72 69 6e 74 66 27 | 73 20 74 68 61 74 20 69 | printf'|s that i|
|00000480| 6e 74 3d 3d 73 68 6f 72 | 74 2e 0d 2a 2f 0d 23 69 |nt==shor|t..*/.#i|
|00000490| 6e 63 6c 75 64 65 20 22 | 56 69 64 65 6f 54 6f 6f |nclude "|VideoToo|
|000004a0| 6c 62 6f 78 2e 68 22 0d | 23 64 65 66 69 6e 65 20 |lbox.h".|#define |
|000004b0| 51 55 45 53 54 20 30 0d | 23 69 66 20 51 55 45 53 |QUEST 0.|#if QUES|
|000004c0| 54 0d 09 23 69 6e 63 6c | 75 64 65 20 22 51 75 65 |T..#incl|ude "Que|
|000004d0| 73 74 2e 68 22 0d 23 65 | 6e 64 69 66 0d 23 69 6e |st.h".#e|ndif.#in|
|000004e0| 63 6c 75 64 65 20 22 51 | 75 69 63 6b 33 2e 68 22 |clude "Q|uick3.h"|
|000004f0| 0d 23 69 6e 63 6c 75 64 | 65 20 3c 74 69 6d 65 2e |.#includ|e <time.|
|00000500| 68 3e 0d 23 69 6e 63 6c | 75 64 65 20 22 6e 72 2e |h>.#incl|ude "nr.|
|00000510| 68 22 09 09 2f 2a 20 4e | 75 6d 65 72 69 63 61 6c |h"../* N|umerical|
|00000520| 20 52 65 63 69 70 65 73 | 20 69 6e 20 43 20 2a 2f | Recipes| in C */|
|00000530| 0d 0d 76 6f 69 64 20 6d | 61 69 6e 28 76 6f 69 64 |..void m|ain(void|
|00000540| 29 0d 7b 0d 09 64 61 74 | 61 52 65 63 6f 72 64 20 |).{..dat|aRecord |
|00000550| 64 61 74 61 2c 6d 6f 6e | 6f 74 6f 6e 69 63 44 61 |data,mon|otonicDa|
|00000560| 74 61 3b 0d 09 63 6f 6e | 74 72 61 73 74 52 65 63 |ta;..con|trastRec|
|00000570| 6f 72 64 20 2a 63 50 74 | 72 3b 0d 09 70 61 72 61 |ord *cPt|r;..para|
|00000580| 6d 52 65 63 6f 72 64 20 | 70 61 72 61 6d 73 3b 0d |mRecord |params;.|
|00000590| 09 69 6e 74 20 69 3b 0d | 09 64 6f 75 62 6c 65 20 |.int i;.|.double |
|000005a0| 63 68 69 53 71 75 61 72 | 65 2c 77 65 69 62 75 6c |chiSquar|e,weibul|
|000005b0| 6c 4c 4c 2c 6d 6f 6e 6f | 74 6f 6e 69 63 4c 4c 3b |lLL,mono|tonicLL;|
|000005c0| 09 2f 2a 20 6c 6f 67 20 | 6c 69 6b 65 6c 69 68 6f |./* log |likeliho|
|000005d0| 6f 64 20 2a 2f 0d 09 69 | 6e 74 20 63 68 69 53 71 |od */..i|nt chiSq|
|000005e0| 75 61 72 65 44 46 2c 77 | 65 69 62 75 6c 6c 44 46 |uareDF,w|eibullDF|
|000005f0| 2c 6d 6f 6e 6f 74 6f 6e | 69 63 44 46 3b 09 2f 2a |,monoton|icDF;./*|
|00000600| 20 64 65 67 72 65 65 73 | 20 6f 66 20 66 72 65 65 | degrees| of free|
|00000610| 64 6f 6d 20 2a 2f 0d 09 | 64 6f 75 62 6c 65 20 70 |dom */..|double p|
|00000620| 2c 6c 6f 67 43 2c 72 61 | 6e 67 65 3b 0d 09 6c 6f |,logC,ra|nge;..lo|
|00000630| 6e 67 20 69 64 75 6d 3b | 0d 09 69 6e 74 20 63 6f |ng idum;|..int co|
|00000640| 6e 64 3d 30 3b 0d 09 23 | 69 66 20 51 55 45 53 54 |nd=0;..#|if QUEST|
|00000650| 0d 09 09 51 75 65 73 74 | 20 2a 71 3b 0d 09 23 65 |...Quest| *q;..#e|
|00000660| 6e 64 69 66 0d 09 0d 09 | 52 65 71 75 69 72 65 28 |ndif....|Require(|
|00000670| 30 29 3b 0d 09 70 61 72 | 61 6d 73 2e 6c 6f 67 41 |0);..par|ams.logA|
|00000680| 6c 70 68 61 3d 30 3b 0d | 09 70 61 72 61 6d 73 2e |lpha=0;.|.params.|
|00000690| 62 65 74 61 3d 33 3b 0d | 09 70 61 72 61 6d 73 2e |beta=3;.|.params.|
|000006a0| 67 61 6d 6d 61 3d 30 2e | 31 31 31 31 37 33 3b 20 |gamma=0.|111173; |
|000006b0| 0d 09 70 61 72 61 6d 73 | 2e 64 65 6c 74 61 3d 30 |..params|.delta=0|
|000006c0| 2e 30 31 3b 0d 09 0d 09 | 69 64 75 6d 3d 74 69 6d |.01;....|idum=tim|
|000006d0| 65 28 4e 55 4c 4c 29 3b | 0d 09 70 72 69 6e 74 66 |e(NULL);|..printf|
|000006e0| 28 22 52 61 6e 64 6f 6d | 20 73 65 65 64 20 25 6c |("Random| seed %l|
|000006f0| 64 5c 6e 22 2c 69 64 75 | 6d 29 3b 09 2f 2a 20 53 |d\n",idu|m);./* S|
|00000700| 6f 20 77 65 20 63 61 6e | 20 72 65 70 72 6f 64 75 |o we can| reprodu|
|00000710| 63 65 20 69 6e 74 65 72 | 65 73 74 69 6e 67 20 63 |ce inter|esting c|
|00000720| 61 73 65 73 20 2a 2f 0d | 09 64 61 74 61 2e 63 6f |ases */.|.data.co|
|00000730| 6e 74 72 61 73 74 73 3d | 35 3b 0d 09 72 61 6e 67 |ntrasts=|5;..rang|
|00000740| 65 3d 30 2e 35 3b 09 09 | 09 09 09 09 09 2f 2a 20 |e=0.5;..|...../* |
|00000750| 6c 6f 67 20 63 6f 6e 74 | 72 61 73 74 20 72 61 6e |log cont|rast ran|
|00000760| 67 65 2c 20 63 65 6e 74 | 65 72 65 64 20 6f 6e 20 |ge, cent|ered on |
|00000770| 6c 6f 67 41 6c 70 68 61 | 20 2a 2f 0d 09 72 61 6e |logAlpha| */..ran|
|00000780| 67 65 3d 32 3b 0d 09 66 | 6f 72 28 69 3d 30 3b 69 |ge=2;..f|or(i=0;i|
|00000790| 3c 64 61 74 61 2e 63 6f | 6e 74 72 61 73 74 73 3b |<data.co|ntrasts;|
|000007a0| 69 2b 2b 29 7b 0d 09 09 | 6c 6f 67 43 3d 70 61 72 |i++){...|logC=par|
|000007b0| 61 6d 73 2e 6c 6f 67 41 | 6c 70 68 61 20 2b 20 72 |ams.logA|lpha + r|
|000007c0| 61 6e 67 65 2a 28 69 2f | 28 64 61 74 61 2e 63 6f |ange*(i/|(data.co|
|000007d0| 6e 74 72 61 73 74 73 2d | 31 2e 30 29 20 2d 20 30 |ntrasts-|1.0) - 0|
|000007e0| 2e 35 29 3b 0d 09 09 64 | 61 74 61 2e 63 5b 69 5d |.5);...d|ata.c[i]|
|000007f0| 2e 63 6f 6e 74 72 61 73 | 74 3d 70 6f 77 28 31 30 |.contras|t=pow(10|
|00000800| 2e 30 2c 6c 6f 67 43 29 | 3b 0d 09 09 64 61 74 61 |.0,logC)|;...data|
|00000810| 2e 63 5b 69 5d 2e 74 72 | 69 61 6c 73 3d 31 30 3b |.c[i].tr|ials=10;|
|00000820| 09 09 09 09 2f 2a 20 74 | 72 69 61 6c 73 20 61 74 |..../* t|rials at|
|00000830| 20 74 68 69 73 20 63 6f | 6e 74 72 61 73 74 20 2a | this co|ntrast *|
|00000840| 2f 0d 09 09 64 61 74 61 | 2e 63 5b 69 5d 2e 63 6f |/...data|.c[i].co|
|00000850| 72 72 65 63 74 3d 62 6e | 6c 64 65 76 28 57 65 69 |rrect=bn|ldev(Wei|
|00000860| 62 75 6c 6c 28 64 61 74 | 61 2e 63 5b 69 5d 2e 63 |bull(dat|a.c[i].c|
|00000870| 6f 6e 74 72 61 73 74 2c | 26 70 61 72 61 6d 73 29 |ontrast,|&params)|
|00000880| 2c 64 61 74 61 2e 63 5b | 69 5d 2e 74 72 69 61 6c |,data.c[|i].trial|
|00000890| 73 2c 26 69 64 75 6d 29 | 3b 0d 09 7d 0d 09 23 69 |s,&idum)|;..}..#i|
|000008a0| 66 20 51 55 45 53 54 0d | 09 09 2f 2a 20 51 75 65 |f QUEST.|../* Que|
|000008b0| 73 74 20 50 61 72 61 6d | 65 74 65 72 73 20 2a 2f |st Param|eters */|
|000008c0| 0d 09 09 71 3d 28 51 75 | 65 73 74 20 2a 29 6d 61 |...q=(Qu|est *)ma|
|000008d0| 6c 6c 6f 63 28 73 69 7a | 65 6f 66 28 51 75 65 73 |lloc(siz|eof(Ques|
|000008e0| 74 29 29 3b 0d 09 09 61 | 73 73 65 72 74 28 71 21 |t));...a|ssert(q!|
|000008f0| 3d 4e 55 4c 4c 29 3b 0d | 09 09 71 2d 3e 6e 43 6f |=NULL);.|..q->nCo|
|00000900| 6e 64 73 3d 31 3b 09 0d | 09 09 71 2d 3e 6e 4c 65 |nds=1;..|..q->nLe|
|00000910| 76 65 6c 73 3d 36 30 30 | 3b 0d 09 09 71 2d 3e 6e |vels=600|;...q->n|
|00000920| 54 72 69 61 6c 73 3d 30 | 3b 0d 09 09 71 2d 3e 67 |Trials=0|;...q->g|
|00000930| 72 61 69 6e 3d 30 2e 30 | 31 3b 09 2f 2a 20 73 74 |rain=0.0|1;./* st|
|00000940| 65 70 20 73 69 7a 65 20 | 6f 66 20 67 72 69 64 2c |ep size |of grid,|
|00000950| 20 69 6e 20 6c 6f 67 20 | 63 6f 6e 74 72 61 73 74 | in log |contrast|
|00000960| 20 2a 2f 0d 09 09 71 2d | 3e 69 6e 69 74 69 61 6c | */...q-|>initial|
|00000970| 53 44 3d 31 3b 0d 09 20 | 09 71 2d 3e 6e 52 65 73 |SD=1;.. |.q->nRes|
|00000980| 70 6f 6e 73 65 73 3d 32 | 3b 0d 09 20 09 71 2d 3e |ponses=2|;.. .q->|
|00000990| 71 75 61 6e 74 69 6c 65 | 4f 72 64 65 72 3d 4e 41 |quantile|Order=NA|
|000009a0| 4e 3b 0d 09 20 09 71 2d | 3e 66 61 6b 65 49 74 3d |N;.. .q-|>fakeIt=|
|000009b0| 30 3b 0d 09 20 09 71 2d | 3e 66 75 6e 63 74 69 6f |0;.. .q-|>functio|
|000009c0| 6e 3d 57 65 69 62 75 6c | 6c 50 52 65 73 70 6f 6e |n=Weibul|lPRespon|
|000009d0| 73 65 3b 0d 09 09 71 2d | 3e 62 65 74 61 3d 70 61 |se;...q-|>beta=pa|
|000009e0| 72 61 6d 73 2e 62 65 74 | 61 3b 0d 09 09 71 2d 3e |rams.bet|a;...q->|
|000009f0| 67 61 6d 6d 61 3d 70 61 | 72 61 6d 73 2e 67 61 6d |gamma=pa|rams.gam|
|00000a00| 6d 61 3b 0d 09 09 71 2d | 3e 64 65 6c 74 61 3d 70 |ma;...q-|>delta=p|
|00000a10| 61 72 61 6d 73 2e 64 65 | 6c 74 61 3b 0d 09 09 71 |arams.de|lta;...q|
|00000a20| 2d 3e 65 70 73 69 6c 6f | 6e 3d 30 2e 30 3b 0d 09 |->epsilo|n=0.0;..|
|00000a30| 09 63 6f 6e 64 3d 30 3b | 0d 09 09 71 2d 3e 67 75 |.cond=0;|...q->gu|
|00000a40| 65 73 73 5b 63 6f 6e 64 | 5d 3d 30 3b 0d 09 09 51 |ess[cond|]=0;...Q|
|00000a50| 75 65 73 74 4f 70 65 6e | 28 71 29 3b 0d 09 09 66 |uestOpen|(q);...f|
|00000a60| 6f 72 28 69 3d 30 3b 69 | 3c 64 61 74 61 2e 63 6f |or(i=0;i|<data.co|
|00000a70| 6e 74 72 61 73 74 73 3b | 69 2b 2b 29 7b 0d 09 09 |ntrasts;|i++){...|
|00000a80| 09 6c 6f 67 43 3d 6c 6f | 67 31 30 28 64 61 74 61 |.logC=lo|g10(data|
|00000a90| 2e 63 5b 69 5d 2e 63 6f | 6e 74 72 61 73 74 29 3b |.c[i].co|ntrast);|
|00000aa0| 0d 09 09 09 66 6f 72 28 | 6a 3d 30 3b 6a 3c 64 61 |....for(|j=0;j<da|
|00000ab0| 74 61 2e 63 5b 69 5d 2e | 63 6f 72 72 65 63 74 3b |ta.c[i].|correct;|
|00000ac0| 6a 2b 2b 29 51 75 65 73 | 74 55 70 64 61 74 65 28 |j++)Ques|tUpdate(|
|00000ad0| 71 2c 63 6f 6e 64 2c 6c | 6f 67 43 2c 31 29 3b 0d |q,cond,l|ogC,1);.|
|00000ae0| 09 09 09 66 6f 72 28 3b | 6a 3c 64 61 74 61 2e 63 |...for(;|j<data.c|
|00000af0| 5b 69 5d 2e 74 72 69 61 | 6c 73 3b 6a 2b 2b 29 51 |[i].tria|ls;j++)Q|
|00000b00| 75 65 73 74 55 70 64 61 | 74 65 28 71 2c 63 6f 6e |uestUpda|te(q,con|
|00000b10| 64 2c 6c 6f 67 43 2c 30 | 29 3b 0d 09 09 7d 0d 09 |d,logC,0|);...}..|
|00000b20| 09 71 5f 72 65 6d 6f 76 | 65 50 72 69 6f 72 28 71 |.q_remov|ePrior(q|
|00000b30| 2d 3e 71 43 6f 6e 64 73 | 5b 63 6f 6e 64 5d 29 3b |->qConds|[cond]);|
|00000b40| 0d 09 09 6d 6f 64 65 20 | 3d 20 71 5f 6d 6f 64 65 |...mode |= q_mode|
|00000b50| 28 71 2d 3e 71 43 6f 6e | 64 73 5b 63 6f 6e 64 5d |(q->qCon|ds[cond]|
|00000b60| 29 3b 0d 09 09 51 75 65 | 73 74 43 6c 6f 73 65 28 |);...Que|stClose(|
|00000b70| 71 29 3b 0d 09 23 65 6e | 64 69 66 0d 09 70 72 69 |q);..#en|dif..pri|
|00000b80| 6e 74 66 28 22 54 65 73 | 74 69 6e 67 20 74 68 65 |ntf("Tes|ting the|
|00000b90| 20 66 75 6e 63 74 69 6f | 6e 20 50 73 79 63 68 6f | functio|n Psycho|
|00000ba0| 6d 65 74 72 69 63 46 69 | 74 2e 5c 6e 22 29 3b 0d |metricFi|t.\n");.|
|00000bb0| 09 70 72 69 6e 74 66 28 | 22 53 69 6d 75 6c 61 74 |.printf(|"Simulat|
|00000bc0| 69 6e 67 20 61 6e 20 6f | 62 73 65 72 76 65 72 20 |ing an o|bserver |
|00000bd0| 77 69 74 68 20 61 20 57 | 65 69 62 75 6c 6c 20 70 |with a W|eibull p|
|00000be0| 73 79 63 68 6f 6d 65 74 | 72 69 63 20 66 75 6e 63 |sychomet|ric func|
|00000bf0| 74 69 6f 6e 2e 5c 6e 22 | 29 3b 0d 09 77 65 69 62 |tion.\n"|);..weib|
|00000c00| 75 6c 6c 44 46 3d 31 3b | 09 2f 2a 20 6e 75 6d 62 |ullDF=1;|./* numb|
|00000c10| 65 72 20 6f 66 20 70 61 | 72 61 6d 65 74 65 72 73 |er of pa|rameters|
|00000c20| 20 74 6f 20 62 65 20 61 | 64 6a 75 73 74 65 64 20 | to be a|djusted |
|00000c30| 69 6e 20 66 69 74 74 69 | 6e 67 20 2a 2f 0d 09 70 |in fitti|ng */..p|
|00000c40| 72 69 6e 74 66 28 22 54 | 68 65 20 73 69 6d 75 6c |rintf("T|he simul|
|00000c50| 61 74 65 64 20 64 61 74 | 61 20 77 69 6c 6c 20 62 |ated dat|a will b|
|00000c60| 65 20 66 69 74 20 75 73 | 69 6e 67 20 25 64 20 64 |e fit us|ing %d d|
|00000c70| 65 67 72 65 65 73 20 6f | 66 20 66 72 65 65 64 6f |egrees o|f freedo|
|00000c80| 6d 2e 5c 6e 22 2c 28 69 | 6e 74 29 77 65 69 62 75 |m.\n",(i|nt)weibu|
|00000c90| 6c 6c 44 46 29 3b 0d 09 | 70 72 69 6e 74 66 28 22 |llDF);..|printf("|
|00000ca0| 4f 62 73 65 72 76 65 72 | 3a 20 22 29 3b 0d 09 70 |Observer|: ");..p|
|00000cb0| 72 69 6e 74 66 28 22 6c | 6f 67 41 6c 70 68 61 25 |rintf("l|ogAlpha%|
|00000cc0| 36 2e 32 66 2c 20 62 65 | 74 61 25 34 2e 31 66 2c |6.2f, be|ta%4.1f,|
|00000cd0| 20 67 61 6d 6d 61 25 35 | 2e 32 66 2c 20 64 65 6c | gamma%5|.2f, del|
|00000ce0| 74 61 25 35 2e 32 66 5c | 6e 22 2c 70 61 72 61 6d |ta%5.2f\|n",param|
|00000cf0| 73 2e 6c 6f 67 41 6c 70 | 68 61 2c 70 61 72 61 6d |s.logAlp|ha,param|
|00000d00| 73 2e 62 65 74 61 2c 70 | 61 72 61 6d 73 2e 67 61 |s.beta,p|arams.ga|
|00000d10| 6d 6d 61 2c 70 61 72 61 | 6d 73 2e 64 65 6c 74 61 |mma,para|ms.delta|
|00000d20| 29 3b 0d 09 70 3d 50 73 | 79 63 68 6f 6d 65 74 72 |);..p=Ps|ychometr|
|00000d30| 69 63 46 69 74 28 26 70 | 61 72 61 6d 73 2c 26 57 |icFit(&p|arams,&W|
|00000d40| 65 69 62 75 6c 6c 2c 26 | 64 61 74 61 2c 26 77 65 |eibull,&|data,&we|
|00000d50| 69 62 75 6c 6c 4c 4c 2c | 77 65 69 62 75 6c 6c 44 |ibullLL,|weibullD|
|00000d60| 46 2c 26 63 68 69 53 71 | 75 61 72 65 2c 26 63 68 |F,&chiSq|uare,&ch|
|00000d70| 69 53 71 75 61 72 65 44 | 46 29 3b 0d 09 70 72 69 |iSquareD|F);..pri|
|00000d80| 6e 74 66 28 22 46 69 74 | 3a 20 20 20 20 20 20 22 |ntf("Fit|:      "|
|00000d90| 29 3b 0d 09 70 72 69 6e | 74 66 28 22 6c 6f 67 41 |);..prin|tf("logA|
|00000da0| 6c 70 68 61 25 36 2e 32 | 66 2c 20 62 65 74 61 25 |lpha%6.2|f, beta%|
|00000db0| 34 2e 31 66 2c 20 67 61 | 6d 6d 61 25 35 2e 32 66 |4.1f, ga|mma%5.2f|
|00000dc0| 2c 20 64 65 6c 74 61 25 | 35 2e 32 66 2c 20 73 69 |, delta%|5.2f, si|
|00000dd0| 67 6e 69 66 69 63 61 6e | 63 65 25 35 2e 32 66 5c |gnifican|ce%5.2f\|
|00000de0| 6e 22 2c 70 61 72 61 6d | 73 2e 6c 6f 67 41 6c 70 |n",param|s.logAlp|
|00000df0| 68 61 2c 70 61 72 61 6d | 73 2e 62 65 74 61 2c 70 |ha,param|s.beta,p|
|00000e00| 61 72 61 6d 73 2e 67 61 | 6d 6d 61 2c 70 61 72 61 |arams.ga|mma,para|
|00000e10| 6d 73 2e 64 65 6c 74 61 | 2c 70 29 3b 0d 09 23 69 |ms.delta|,p);..#i|
|00000e20| 66 20 51 55 45 53 54 0d | 09 09 70 72 69 6e 74 66 |f QUEST.|..printf|
|00000e30| 28 22 51 55 45 53 54 20 | 6d 6f 64 65 20 25 2e 32 |("QUEST |mode %.2|
|00000e40| 66 5c 6e 22 2c 6d 6f 64 | 65 29 3b 0d 09 23 65 6e |f\n",mod|e);..#en|
|00000e50| 64 69 66 0d 0d 09 2f 2a | 0d 09 57 65 27 72 65 20 |dif.../*|..We're |
|00000e60| 64 6f 6e 65 2c 20 62 75 | 74 20 6a 75 73 74 20 74 |done, bu|t just t|
|00000e70| 6f 20 73 68 6f 77 20 6f | 66 66 2c 20 6c 65 74 27 |o show o|ff, let'|
|00000e80| 73 20 70 72 69 6e 74 20 | 6f 75 74 20 65 76 65 72 |s print |out ever|
|00000e90| 79 74 68 69 6e 67 20 74 | 68 61 74 20 61 6e 79 6f |ything t|hat anyo|
|00000ea0| 6e 65 0d 09 63 6f 75 6c | 64 20 70 6f 73 73 69 62 |ne..coul|d possib|
|00000eb0| 6c 79 20 77 61 6e 74 2e | 20 49 6e 20 72 65 61 6c |ly want.| In real|
|00000ec0| 20 6c 69 66 65 20 49 20 | 77 6f 75 6c 64 20 73 6b | life I |would sk|
|00000ed0| 69 70 20 74 68 69 73 20 | 6a 75 6e 6b 2e 0d 09 2a |ip this |junk...*|
|00000ee0| 2f 0d 09 6d 6f 6e 6f 74 | 6f 6e 69 63 44 61 74 61 |/..monot|onicData|
|00000ef0| 3d 64 61 74 61 3b 0d 09 | 4d 6f 6e 6f 74 6f 6e 69 |=data;..|Monotoni|
|00000f00| 63 46 69 74 28 26 6d 6f | 6e 6f 74 6f 6e 69 63 44 |cFit(&mo|notonicD|
|00000f10| 61 74 61 2c 26 6d 6f 6e | 6f 74 6f 6e 69 63 4c 4c |ata,&mon|otonicLL|
|00000f20| 2c 26 6d 6f 6e 6f 74 6f | 6e 69 63 44 46 29 3b 09 |,&monoto|nicDF);.|
|00000f30| 2f 2a 20 6f 76 65 72 77 | 72 69 74 65 73 20 64 61 |/* overw|rites da|
|00000f40| 74 61 20 77 69 74 68 20 | 66 69 74 20 2a 2f 0d 09 |ta with |fit */..|
|00000f50| 70 72 69 6e 74 66 28 22 | 5c 6e 63 6f 6e 74 72 61 |printf("|\ncontra|
|00000f60| 73 74 20 54 72 69 61 6c | 73 20 52 69 67 68 74 20 |st Trial|s Right |
|00000f70| 20 52 61 74 69 6f 20 20 | 57 65 69 62 75 6c 6c 20 | Ratio  |Weibull |
|00000f80| 4d 6f 6e 6f 74 6f 6e 65 | 5c 6e 22 29 3b 0d 09 66 |Monotone|\n");..f|
|00000f90| 6f 72 28 69 3d 30 3b 69 | 3c 64 61 74 61 2e 63 6f |or(i=0;i|<data.co|
|00000fa0| 6e 74 72 61 73 74 73 3b | 69 2b 2b 29 7b 0d 09 09 |ntrasts;|i++){...|
|00000fb0| 63 50 74 72 3d 26 64 61 | 74 61 2e 63 5b 69 5d 3b |cPtr=&da|ta.c[i];|
|00000fc0| 0d 09 09 70 72 69 6e 74 | 66 28 22 25 36 2e 33 66 |...print|f("%6.3f|
|00000fd0| 20 20 20 25 35 6c 64 20 | 25 35 6c 64 20 25 37 2e |   %5ld |%5ld %7.|
|00000fe0| 33 66 20 25 37 2e 33 66 | 20 25 37 2e 33 66 5c 6e |3f %7.3f| %7.3f\n|
|00000ff0| 22 2c 0d 09 09 09 63 50 | 74 72 2d 3e 63 6f 6e 74 |",....cP|tr->cont|
|00001000| 72 61 73 74 2c 63 50 74 | 72 2d 3e 74 72 69 61 6c |rast,cPt|r->trial|
|00001010| 73 2c 63 50 74 72 2d 3e | 63 6f 72 72 65 63 74 2c |s,cPtr->|correct,|
|00001020| 0d 09 09 09 63 50 74 72 | 2d 3e 63 6f 72 72 65 63 |....cPtr|->correc|
|00001030| 74 2f 28 64 6f 75 62 6c | 65 29 63 50 74 72 2d 3e |t/(doubl|e)cPtr->|
|00001040| 74 72 69 61 6c 73 2c 0d | 09 09 09 57 65 69 62 75 |trials,.|...Weibu|
|00001050| 6c 6c 28 63 50 74 72 2d | 3e 63 6f 6e 74 72 61 73 |ll(cPtr-|>contras|
|00001060| 74 2c 26 70 61 72 61 6d | 73 29 2c 0d 09 09 09 6d |t,&param|s),....m|
|00001070| 6f 6e 6f 74 6f 6e 69 63 | 44 61 74 61 2e 63 5b 69 |onotonic|Data.c[i|
|00001080| 5d 2e 63 6f 72 72 65 63 | 74 2f 28 64 6f 75 62 6c |].correc|t/(doubl|
|00001090| 65 29 6d 6f 6e 6f 74 6f | 6e 69 63 44 61 74 61 2e |e)monoto|nicData.|
|000010a0| 63 5b 69 5d 2e 74 72 69 | 61 6c 73 0d 09 09 09 29 |c[i].tri|als....)|
|000010b0| 3b 0d 09 7d 0d 0d 09 63 | 68 69 53 71 75 61 72 65 |;..}...c|hiSquare|
|000010c0| 3d 32 2e 30 2a 28 6d 6f | 6e 6f 74 6f 6e 69 63 4c |=2.0*(mo|notonicL|
|000010d0| 4c 2d 77 65 69 62 75 6c | 6c 4c 4c 29 3b 0d 09 63 |L-weibul|lLL);..c|
|000010e0| 68 69 53 71 75 61 72 65 | 44 46 3d 6d 6f 6e 6f 74 |hiSquare|DF=monot|
|000010f0| 6f 6e 69 63 44 46 2d 77 | 65 69 62 75 6c 6c 44 46 |onicDF-w|eibullDF|
|00001100| 3b 0d 09 70 3d 50 43 68 | 69 53 71 75 61 72 65 28 |;..p=PCh|iSquare(|
|00001110| 63 68 69 53 71 75 61 72 | 65 2c 63 68 69 53 71 75 |chiSquar|e,chiSqu|
|00001120| 61 72 65 44 46 29 3b 0d | 09 70 72 69 6e 74 66 28 |areDF);.|.printf(|
|00001130| 22 5c 6e 43 68 69 20 73 | 71 75 61 72 65 20 25 2e |"\nChi s|quare %.|
|00001140| 31 66 20 77 69 74 68 20 | 25 64 20 64 65 67 72 65 |1f with |%d degre|
|00001150| 65 73 20 6f 66 20 66 72 | 65 65 64 6f 6d 2c 20 79 |es of fr|eedom, y|
|00001160| 69 65 6c 64 69 6e 67 20 | 61 20 73 69 67 6e 69 66 |ielding |a signif|
|00001170| 69 63 61 6e 63 65 20 6f | 66 20 25 2e 32 66 5c 6e |icance o|f %.2f\n|
|00001180| 22 0d 09 09 2c 63 68 69 | 53 71 75 61 72 65 2c 28 |"...,chi|Square,(|
|00001190| 69 6e 74 29 63 68 69 53 | 71 75 61 72 65 44 46 2c |int)chiS|quareDF,|
|000011a0| 70 29 3b 0d 7d 00 00 00 | 00 00 00 00 00 00 00 00 |p);.}...|........|
|000011b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001200| 00 00 01 00 00 00 01 54 | 00 00 00 54 00 00 00 46 |.......T|...T...F|
|00001210| 0d 0d 26 8c 10 14 4c ee | 18 00 ff f8 4e 5e 4e 75 |..&...L.|....N^Nu|
|00001220| 87 73 6b 69 70 5f 62 6c | 00 00 4e 56 ff fe 48 e7 |.skip_bl|..NV..H.|
|00001230| 15 54 65 73 74 50 73 79 | 63 68 6f 6d 65 74 72 69 |.TestPsy|chometri|
|00001240| 63 46 69 74 2e 63 02 00 | 00 00 54 45 58 54 43 57 |cFit.c..|..TEXTCW|
|00001250| 00 00 54 45 58 54 43 57 | 49 45 01 00 01 04 00 00 |..TEXTCW|IE......|
|00001260| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001270| 00 00 a2 41 95 96 00 00 | 11 25 00 00 01 9a 0c 47 |...A....|.%.....G|
|00001280| 00 10 6d 02 5d 47 30 06 | c1 c5 d0 47 3c 00 52 94 |..m.]G0.|...G<.R.|
|00001290| 78 01 60 dc 52 47 4a 33 | 70 00 66 d8 4a 04 66 0c |x.`.RGJ3|p.f.J.f.|
|000012a0| 2f 14 48 6d f9 6a 4e ba | f8 10 50 4f 4a 6e ff fe |/.Hm.jN.|..POJn..|
|000012b0| 6c 08 48 c6 20 06 44 80 | 3c 00 48 c6 20 06 4c ee |l.H. .D.|<.H. .L.|
|000012c0| 18 f0 ff e6 4e 5e 4e 75 | 86 73 5f 61 74 6f 69 00 |....N^Nu|.s_atoi.|
|000012d0| 00 00 4e 56 ff f8 48 e7 | 01 18 26 6e 00 0c 28 6e |..NV..H.|..&n..(n|
|000012e0| 00 08 60 2e 10 1c 48 80 | 48 c0 2f 00 4e ad 01 d2 |..`...H.|H./.N...|
|000012f0| 2d 40 ff fc 10 1b 48 80 | 48 c0 2f 00 4e ad 01 d2 |-@....H.|H./.N...|
|00001300| 00 00 00 48 00 09 4d 6f | 6e 61 63 6f 00 00 00 00 |...H..Mo|naco....|
|00001310| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001320| 00 00 00 00 00 00 00 03 | 00 04 00 3d 00 04 01 9c |........|...=....|
|00001330| 02 7c 00 3d 00 04 01 9c | 02 7c ac 07 e1 03 00 00 |.|.=....|.|......|
|00001340| 04 71 00 00 04 71 00 00 | 01 fd 01 00 00 00 00 04 |.q...q..|........|
|00001350| 00 01 00 01 00 00 01 00 | 00 00 01 54 00 00 00 54 |........|...T...T|
|00001360| 00 00 00 46 01 9f 30 88 | 20 b0 00 00 00 1c 00 46 |...F..0.| ......F|
|00001370| 00 01 4d 50 53 52 00 00 | 00 12 4d 57 42 42 00 00 |..MPSR..|..MWBB..|
|00001380| 00 1e 03 ed ff ff 00 00 | 00 00 00 00 00 00 03 f0 |........|........|
|00001390| ff ff 00 00 00 4c 00 00 | 00 00 00 00 00 00 00 00 |.....L..|........|
|000013a0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000013b0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000013c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000013d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000013e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000013f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
+--------+-------------------------+-------------------------+--------+--------+